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ABSTRACT 


The two-layer anisotropic random medium model is used to investigate fully polan- 
metric scattering properties of earth terrain media. The polarization covariance matrices 
for the untilted and tilted uniaxial random medium are evaluated using the strong fluc- 
tuation theory and distorted Bom approximation. In order to account for the azimuthal 
randomness in the growth direction of leaves in tree and grass fields, an averaging scheme 
over the azimuthal direction is also applied. It is found that characteristics of terrain 
clutter can be identified through the analysis of each element of the covariance matrix. 
Theoretical results are illustrated by the comparison with experimental data provided by 
MIT Lincoln Laboratory for tree and grass fields. 
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I. INTRODUCTION 


The layered random medium model for terrain cover, in conjunction with the ap- 
plication of electromagnetic wave theory, provides a systematic approach in relating the 
radar backscatter response to physical properties of geophysical media, where the volume 
scattering and anisotropic effects are attributed to the embedded inhomogeneities with 
elongated geometry and preferred alignment. For active and passive microwave remote 
sensing, the effectiveness of this model has been demonstrated by extraction of physical 
parameters from data matching for scene interpretation and feature identification [Zuniga 
et al., 1980; Tsang and Kong, 1981; Tsang et al., 1982; Lee and Kong, 1985; Lin et al., 
1988a, b]. In polarimetric microwave remote sensing, the random medium model has also 
been applied to radar image simulation, terrain clutter identification and classification, 
and radar range profile simulation of terrain clutter [Shin et al., 1986; Borgeaud et al., 
1987a; Borgeaud, 1987b; Kong et al., 1988; Yueh et al., 1988a; Swartz et al . , 1988]. 

The theoretical modeling of earth terrain media is described in Section II. To account 
for the large permittivity contrast between the background medium and the embedded in- 
homogeneities, the strong fluctuation theory is applied in Section III, in order to derive 
the external field which is decomposed into mean (coherent part) and scattered field (in- 
coherent part) components. The mean field is derived using the Feynman diagrammatic 
technique and bilocal approximation. In Section IV, the effective permittivity tensor of 
the medium is obtained from the constitutive relation. The effective permittivity tensor 
is computed in the low frequency limit for an unbounded uniaxial random medium with 
two-phase mixture. The distorted Born approximation is then used to calculate the scat- 
tered field in Section V. The polarization covariance matrices for the untilted and tilted 
uniaxial random medium, as well as that obtained from the azimuthal averaging of the 
tilted case are discussed in Section VI. Theoretical' calculations of the fully polarimetric 
random medium model are compared with the experimental data for trees and a grass field 
supplied by MIT Lincoln Laboratory. 

II. TWO-LAYER ANISOTROPIC RANDOM MEDIUM MODEL 

The two-layer ani sotropic random medium model is applied to simulate earth terrain 
media such as tree canopy and a grass field. In region 1 of Figure 1, the random medium 
with height d is characterized by a spatially random permittivity tensor fi(F). Regions 
0 and 2 (air and ground) are considered to be homogeneous media with permittivities e 0 
and £j, respectively. All three regions are assumed to have the same permeability, /i<>. 
The polarimetric radar system is located in region 0. The coordinate system (z, y, z ) is 



oriented so that the xy plane coinrides with the air/clutter interface. Interfaces at z = 0 
and z = —d are assumed to be planar, extending to infinity on the horizontal plane, and 
parallel to each other. The volume scattering effect is caused by randomly distributed 
scatterers embedded in the layer (e.g., moisture content in leaves) while the anisotropic 
effect is due to the elongated shape and preferred alignment of the scatterers. As shown in 
Figure 2, the moisture content in trees and grass is modeled as cylindrical scatterers with 
the preferred alignment in the yz plane. The elongated direction is along the z'-axis and 
tilted by an angle ip with respect to the z-axis. Hence, the coordinate system (x', y', z') 
which describes the orientation of the scatterers is rotated clockwise by an angle ip about 
the z-axis. In this manner, the z'-axis can be chosen as the optic axis of the permittivity 
tensor (r) [Kong, 1986]. 

In the polarimetric microwave remote sensing of terrain cover, the polarization 
covariance matrix is essential for terrain-feature identification and classification [Borgeaud 
ct al., 1987a; Kong et ai, 1988]. For a plane wave impinging on the random medium layer, 
the scattered electric field E 0 (f) is related to the incident electric field ~E 0 j(F) by 
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where the horizontal and vertical components of the the incident and scattered electric 
fields are expressed as £\,, E Kl , E, t , and 2?„, respectively. For reciprocal media, the 
relation, u = /*», is obtained in the backscattering direction. Thus, the covariance 
matrix is defined as 
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where A is the illuminated area and denotes the ensemble average of the random 

variable X. The covariance matrix is normalized so that the diagonal terms represent 
the backscattering cross sections per unit area for the HH, HV, and VV polarizations, 
respectively. 
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m. STRONG FLUCTUATION THEORY 

Consider a linearly polarized time-harmonic electromagnetic plane wave impinging 
on the random medium. The time-harmonic factor e“*"‘ is used. The propagation vector 
Km at the observation point has the incident angle So, with respect to the z-axis (Fig. 

1) and the azimuthal angle • The total electric fields in regions 0 and 1 satisfy the 

following vector wave equations: 

V x V x E 0 (r) - A£E 0 (r) = 0 (3) 

V x V x Z t (F) - kf-lQ . = 0 (4) 

where = u> J ^>€o- In the microwave frequency range, the dielectric constants of air and 
moisture content are very distinct. In order to account for the large permittivity contrast 
between the host medium and embedded inhomogeneities, a deterministic quantity ? (1 is 
introduced in the strong fluctuation theory [Ryzhov and Tamoikin, 1970; Tseng and Kong, 
1981; Tsang et al., 1982; Lin et al., 1988a, b]. Equation (4) can then be rewritten as 

V x V x E.(?) ®i(f) - t'o ( ?l(f) < ~ ? " ) • ®i00 

= • £,(F) (5) 

Physically, f #1 is the effective permittivity tensor of the medium in the low frequency 
range where the scattering loss is negligible because the size of scatterers is much smaller 
than the wavelength of the incident field [Tsang et al., 1985]. Treating the term on the 
right-hand side of (5) as the effective source, the total electric fields in regions 0 and 1 can 
be represented in integral forms as follows: 

~E 0 (¥) = ^(r) + *o / ‘friSwOSM * • %i(¥i) (6) 

JVt 

Hi(¥) = E ( l ° ) (r) + fc* / d*Fi(?xi(^»^i) • (7) 

where the unperturbed electric fields, and E^ o) (r), are solutions to the homoge- 

neous vector wave equations in (3) and (5) in the absence of the effective source term and 
the subscript Vi denotes the volume of region 1. The dyadic Green’s functions, ^oi( ? > ^i) 
and frn(r, r t ), which correspond to the responses at r in regions 0 and 1, respectively, 
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due to a point source at r x in the anisotropic homogeneous medium with the permittivity 
tensor f #x , are governed by the following vector wave equations: 


V x V x Cx 0 i(r, F x ) - fco^oiO*, r x ) = 0 


V x V x <?„(F, 7,) - • ^ XX (F, 7,) = lS(7, 7,) 

Co 


( 8 ) 


(9) 


When the source and the observation points are in the same region, there is a singularity 
in the dyadic Green’s function, which can be dealt with by separating ^„(7, F x ) into 
principal value PSU„(7, r x ) and singular components, namely, 
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Sn(r, r x ) = PSZ„(7, F x ) - %6(7 - 7,) 

*0 


( 10 ) 


Thus, substituting (10) into (7) yields 

7,(7) = ZJV) + / <P7,PSd„(7 t 7,) • ( X (F X ) • 7,(7,) (11) 

'Vi 

r 1 (F)»[7+S l .9 l (F)].r,(F) ( 12 ) 

fi( ? ) 3 5 i( ? ) • P + 5 t - 5 i(F)] ' ( 13 ) 

Physically, F,(7) can be interpreted as the external field [Ryzhov and Tamoikin, 1970]. 

The external field can be decomposed into a coherent part (mean field) (7,(7)} 
and an incoherent part 7,(7) (scattered field) [Tsang et al., 1985]. Under the bilocal 
approximation, the mean field (7,(7)} can be expressed as 

(7,(7)} = T?, o) (7) + %f Y J Vt d , F x< fF a P5^ xx (F,r x ) -j u , f (7 u 7,) ■ (7,(7,)} (14) 


where 


( 15 ) 


= *.’(«,(?.) • PS77 n (T„T,) ■ {>(*,)) 

Using the index notation, (15) can be expressed as 

[?!.//(*,*)], = klT lmKf R ( (7 u 7,) pn(F 1 ,F 1 )] mw + r, mn ,S^S(7, - 7 ,) (16) 
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where Equation (10) was used and ^£wm(n)£i«p( ? i)y has been replaced with the prod- 
uct of the variance, T tmnr (l,m,n,p = x,y,z), and the normalized correlation function, 
Rf(W lt W 3 ). For a statistically homogeneous medium, the normalized correlation function 
is a function of the displacement between F t and Fj, i.e., 

= ~ ^i) (17) 

Notice that the intrinsic properties of the embedded inhomogeneities are directly related 
to the volume scattering mechanism through the normalized correlation function in (14). 

IV. EFFECTIVE PERMITTIVITY 

The effective permittivity tensor for the uniaxial random medium is derived in this 
section, following the usual approach where the effect of boundary layers is neglected [Tsang 
and Kong, 1981; Stogryn, 1984]. In an unbounded uniaxial medium, with a quasi-static 
effective permittivity tensor of the form 

'€,0 0 

! fl = 0 t, 0 (18) 

.0 0 €„ . 

the mean field satisfies the following vector wave equation 

V x V x (?t(F)) - k* 0 € -£ • ( J\(F )) + * 0 * d , r l f l .„(F,F 1 ) • (7\(F0) = 0 (19) 

The dispersion relation, obtained from (19), allows the effective permittivity tensor for the 
medium to be defined as 

ci 0 O' 

!,.„(*) = = 0 ‘1 0 ( 20 ) 

.0 0 « u . 

where ( uff (k) is the Fourier transform of ? w// (ri - F,) as given in (16). In the low 
frequency limit, £ U f f (k) can be approximated by 

[L„( 0 )]„ = /" [3,(E)]„ + r imnp ^ mn (21) 
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where ♦*(£) i* the Fourier transform of the normalized correlation function The 

Fourier transform of the dyadic Green’s function (7,(r) is given as 

■ k * —k,ky O' 

&,(k) = 7 t -k y k a kl 0 

(*i + - *0%) 

.0 0 0 . 
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in which k* = AsJ + k*. Analytical expressions for £ u /j(0) and the diagonal tensor ^ 
are shown in Appendix A. For a two-phase mixture consisting of a background medium 
and scatterers with permittivities e» and e, , respectively, l fl can be evaluated from the 
following criterion [Tsang and Kong, 1981]: 


(1(f)) = 


From (13) and (23), we obtain 


!» n _ f \ + — ~ — !* — f = q 

Co + 5(c»- e ,) 1 J ‘ )+ eo + Sit.-*,) 1, 


- '•> + - 0 (25) 

where f, is the fractional volume of scatterers. Equations (24) and (25) are the anisotropic 
version of Polder and van Santen’s mixing formula [Polder and van Santen, 1946]. 
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V. DISTORTED BORN APPROXIMATION 

In the strong fluctuation theory, the original random medium is replaced with an 
equivalent random medium which has the permittivity, ? t i. This can be accomplished 
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without altering the characteristic of the volume scattering effect because the effective 
source defined in (5) remains intact, that is, 


5, (r) • E,(f) = *;(,(?)• 7, (F) (26) 

where (12) and (13) have been used. From (6), the total electric field in region 0 can be 
rewritten as 

Zo(r) = Z* *(r) + &o j v <^oi(F, ?i) • £i(Fi) * PiO^i) 

= eJV) + K(f) m 

where Pq(f) is the scattered electric field due to the effective source. By adopting the 
concept of the distorted Born approximation in quantum mechanics [Newton, 1966; Schiff, 
1968), the scatterers are assumed to be embedded in the equivalent medium with the effec- 
tive permittivity, fw//(0) [Lang, 1981], so that the mean field, ^7\(F)^, is used to ap- 
proximate the external field, J\(F), and the mean dyadic Green’s function, ^^oi(F, F t )^, 
is used to replace (7oi(F, r t ). The following vector wave equations: 

V x V x (7\(F)) - • (^(F)^ = 0 (28) 

V x V x (<7 10 (F, F x )) - • (^io(F, FO) = 0 (29) 

in conjunction with the symmetric property of the dyadic Green’s functions, namely, 
^^oi(F, Fj >>-<* xo(^ii ')> , are used to derive the mean field and the mean dyadic 
Green’s function for the scattered electric field which can be written as 

F 0 (f) = kl f <Cf,(3„(f, f,)) • Ilf,) ■ (7,(f,)) (30) 

Physically, the scattered electric field, under the distorted Born approximation, corre- 
sponds to the single scattering of the coherent fields [Tsang, et al. f 1985]. It is also known 
as the first-order multiple scattering [Ishimaru, 1978]. After decomposing ^ 0 (r) hor- 
izontal and vertical components and making use of (1) and (2), the nine elements of the 
polarization covariance matrix can be computed. 


VI. DATA MATCHING AND APPLICATIONS 


Data Matching 


For an untilted uniaxial random medium, the optic axis is in the z-direction, ip = 
Q* For vertically and horizontally polarized incident fields, the single scattering 
process in the distorted Bom approximation does not depolarize the incident field in the 
backscattering direction and the coefficient /*. vanishes. Hence, the covariance matrix 
becomes 


u = 


1 0 Py/j 

0 0 0 


(31) 


o 7 j 

where <r is the backscattering coefficient for the HH polarization, 7 is the copolarization 
intensity ratio (^../cr**), and p is the normalized correlation coefficient between HH and 
VV returns given by 


Oy/7 


( 316 ) 


In general, for the case of tilted uniaxial random medium (if> = jL 0), depo- 
larization effects exist even in the single scattering process, and the nine elements of the 
covariance matrix are all nonzero, i.e., 


where 



Py/l ' 
is/T* 


Lp'y/7 Cyfye 7 J 


(32) 



(32a) 

(32b) 


are the correlation coefficients between the HH and HV and HV and W ch ann els, respec- 
tively. However, when the downward propagation vector *<* lies in the yz plane, i.e., 
<f>oi = 90°, both the double refraction phenomenon and the depolarization effect in the 



backscattering direction disappear so that fk* = 0 [Lee and Kong, 1985a]. Under this 
condition, the covariance matrix also retains the same form as the untilted case, as shown 
in (31). It should be noted that the growth directions of leaves in trees and grass fields are 
not necessarily aligned in one direction. In order to account for the random orientation 
of the scatterers in the azimuthal direction, an azimuthal averaging scheme [Borgeaud, 
1987b] is applied to derive the polarization covariance matrix. It is found that although 
/*• ^ 0, four of the elements in the covariance matrix reduce to zero [Borgeaud et al., 
1987a; Kong et al., 1988]. Thus the polarization covariance matrix becomes 

’ 1 0 Py /i' 

V=<7 0 e 0 (33) 

P*- s/7 0 7 • 

The existence of the four zero elements is due to the fact that the random medium possesses 
an azimuthal symmetry in a statistical sense. The four zero elements of the covariance 
matrix indicate that there is no correlation between the copolarization (HH and W) and 
cross- polarization (HV) returns, i.e., (3 = 0 and ( = 0. 

In order to illustrate the random medium model, polarimetric data for forests and 
grass fields obtained from MIT Lincoln Laboratory are matched with this model. The op- 
erating frequency was 35 GHz, and the angle of incidence was 82* . Measured covariance 
matrices for grass and tree regions, given in Table 1, were calculated from the blocked 
out regions shown in Figure 7 of [Swartz et al., 1988]. Analysis of this database indi- 
cated that there is essentially no correlation between the HH and HV, and between 
the HV and VV polarimetric returns; from a statistical point of view the terrain clutter 
exhibited azimuthal symmetry, and therefore (3 ~ 0 and £ ~ 0 . Thus, the form of the 
covariance matrix shown in (33) was utilized. The corresponding theoretically calculated 
covariance matrix parameters, also shown in Table 1, were obtained by azimuthal averag- 
ing [Borgeaud, 1987b] of the covariance matrices obtained using the two-layer anisotropic 
random medi um model with the application of the strong fluctuation theory and distorted 
Born approximation. 

Applications 

The random medium model provides a systematic approach for the design of optimal 
target detection and classification algorithms [Novak et al., 1987] and the identification of 
terrain media such as vegetation canopy, forest, and snow-covered fields using the optimum 
polarimetric classifier [Kong et al., 1988]. The polarization covariance matrices for various 
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terrain cover* were computed from theoretical model* of random medium. The optimal 
classification scheme made use of a quadratic distance measure and was applied to classify 
a vegetation canopy consisting of both tree* and grass; i.e., using the fully polarimetric 
covariance matrix for the scattered fields, the Baye* likelihood ratio teat was performed on 
specific measurement* to classify the terrain into different categories. The Bayes likelihood 
ratio test has been shown to be optimal in the sense that it minimises the probability of 
error [Fukunaga, 1972; Hord, 1986; Richards, 1986; Swain, 1978]. 

The probabilities of error obtained using fully polarimetric data were compared with 
the probabilities of error obtained using several single, polarimetrically derived features. 
For each single feature studied, the corresponding probability density function was derived. 
Once the probability density functions were known, the Bayes likelihood ratio test was 
performed to classify terrain elements into different categories. For classification schemes 
based on single features, closed-form expressions for the probability of error were calculated 
[Kong et al., 1988]. 

A supervised Bayes classification was applied to synthetic aperture radar (SAR) 
polarimetric images in order to identify their various earth terrain components [Lim et 
al., 1988]. Both fully polarimetric and normalized polarimetric classifications were em- 
ployed to classify radar imagery. It was again shown, in this case throu gh use of radar 
images, that fully polarimetric classifications yielded optimal results; however, an optimal 
normalized classification scheme [Yueh et al., 1988b] indicated improved performance in 
regions where the absolute backscattering coefficients differed significantly from that of the 
training regions due to the variation in return power as a function of incident angle. 

The Bayes classification is known to give minimum probability of error as long as 
the statistical distributions of the radar returns are known and the training areas are 
accurate. However, as noted above, the selected training regions were not sufficient to 
classify correctly the entire image due to the variation in incident angle over the imaged 
swath. By employing the random medium model to generate the classifier training data, 
the effect of the incident angle can be included in the classification scheme. 

Another application of the fully polarimetric, multi-incident angle and frequency 
random medium model is to match training data from various image regions, say at one 
frequency and during one particular season, and then generate classifier training data 
which will predict the backscattered response of the terrain at different frequencies and 
during different seasons of the year. For example, the effects of the winter season can be 
simulated by adding a layer of snow to a terrain from which training data was recorded 
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during the summer season. In this manner, data need not be collected for all seasons and 
at all frequencies. The random medium model can be used to predict the backseat tered 
response based on the initial training data. 


APPENDIX A 

The normalized correlation function, y', z'), for the unbounded uniaxial 

random medium in the (*', y 1 , z') coordinate frame is assumed to be 


exp 
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The components of 3\ are obtained from the frequency-independent part of (A5) and 
(A6): 

5 = 60 I V | W*7» + 1 + q 3 1 

e »(^ + 2)U v /(^ + 1 ) 8 [ v ^I+T-lJ 


^log^ +1 


.V2 




(A8) 


S.= 


_fo_f V 

«#«(*7 a + 2)U, 


2\/2tj ^ f >/2 + 1 


)/(>»' + >) 

v^r]"”’} 


log 


[tjvV + 1 + »? J 


V^TT-l 


(A9) 


*7 = 




N '.v 


(A10) 


From (A8) and (A9), it can be shown that 5 and 5* satisfy the condition [Stogryn, 

1983]: 

2—5 + —S t = 1 (All) 

Co Co 
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Figure 1 : 


Anisotropic two-layer random medium configuration. 



Figure 2 : Geometry of the permittivity tensor in the anisotropic random medium. 
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Table 1 : Covariance matrix elements for grass (a) and trees (b). 
























